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Abstract: We study the matter density fluctuations in the running cosmological constant 
(RCC) model using linear perturbations in the longitudinal gauge. Using this observable 
we calculate the growth rate of structures and the matter power spectrum, and compare 
them with the SDSS data and other available data of the linear growth rate. The dis- 
tribution of collapsed structures may also constraints models of dark energy. It is shown 
that RCC model enhances departures from the ACDM model for both cluster number and 
cumulative cluster number predicted. In general increasing the characteristic parameter 
v leads to significant growth of the cluster number. In general, we found that the theory 
of perturbations provides a good tool to distinguish the new model RCC of the standard 
cosmological model ACDM. 
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1. Introduction 

Recent results from supernova la [[ij and Wilkison Microwave Anisotropy Probe (WMAP), 
especially the five-year data Q indicate that the present Universe is accelerating and that 
has negligible space curvature. 

In this context, the non-relativistic matter contributes about 30% (dark matter plus 
baryonic matter) of the critical density of the Universe and 70% of the energy density 
is not known and is called of dark energy. The dark energy is generally attributed to a 
cosmological constant (CC) and is physically equivalent to vacuum energy. This component 
breaks the strong energy condition, but is the simplest model that one can build. This 
model although satisfactory from an observational point of view, is anyway theoretically 
disfavored, since there is a huge different between the predicted and observed values of 
CC §. Other posibilities have been considered, among these the most popular are based 
on scalar fields, known as quintessence models (4j and models based on scalar fields with 
non-canonical kinetic energy, called of k-essence models Unfortunately, these models 
are difficult of to distinguish from a CC, because the change of the scalar field may be 
extremely slow and there are an degeneration of the potential of quintessence. We may 
also mention other models of dark energy as are: alternative gravity scalar-tensor type 
Q, F(R) theories j?]], and models that introduce extra-dimensions as the DGP model (8|, 
among other. For a recent review of models of dark energy see |J- 
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In practice the dark energy can be seen into of the Friedmann equations through 
its effective energy density and pressure. The ratio of the energy density to pressure, 
known as the equation of state (EOS), ratio w(z) = P(z)/p(z), has important within 
the equations regardless of its physical origin. The function w(z), where z is the redshift, 
is key to understanding the dynamics of the cosmological expansion. Another quantity 
important in the background cosmological is the deceleration parameter q(z). Both are 
the major functions used in studies of constraints of parameters using background data. 
In particular, the deceleration parameter can be used to investigate the kinematics of the 
accelerated expansion using luminosity distances data (and other data, see [1C] and [11]). 

The next step, in the study of a cosmological model, is to use first order cosmological 
perturbations. In fact, the behaviour of linear perturbations in a scalar field and its effect 
on large scale structure formation has been investigated by a number of authors, see [12]. 
Also the behaviour of the nonlinear gravitational collapse has been investigated [13]. These 
studies are fundamentals to understanding and discrimination between competing models. 

In this paper, we investigate the cosmological consequences of an model motivated in 
quantum field theory QFT, more specifically from the renormalization group method. The 
idea to study the renormalization group effects as a way to solve the CC (energy of vacuum) 
has been explored in the papers ||, |l|, II, 0, |l|, ID, H- In particular, a consistent model 
was presented in [^, 22] and is called of Running Cosmological Constant (RCC), where the 
CC is an evolving parameter. In general, in this approach the energy density of the vacuum 
p\(z) can be deduced as a quadratic function of the expansion rate (see equation 2.3 in 
section II), which allows unambiguously resolve the equation of conservation of energy and 
determine the hubble parameter H(z) as function of the redshift. This model is equivalent 
to models of dark energy in the fluid approximation (as are, parameterizations of the EOS 
and quintessence models |23|). The background cosmology of this model has been well 
investigated together with data from supernovae la, restricting the values of the parameter 
v, that represents the "run" of the CC in the RCC model. Also, in [25] was determined 
the matter power spectrum in the synchronous gauge. 

In this paper we focus the analyze of the matter density perturbations using the lon- 
gitudinal gauge, in which Bardeen gauge-invariant variables are identical to the remaining 



metric variables. This approach was pioneered in references [^g, 27]. The formalism is 
applied to our RCC model to determine the linear growth rate and the matter power spec- 
trum as functions of the v parameter. We compared the results with observational data 
from the SDSS (Sloan Digital Sky Survey). 

On the other hand, has long been recognized in various theoretical work [28] and in 
simulations that the evolution of the clusters number (or number counts) can determine 
properties of the dark energy. Therefore, we use our results of linear density perturbations 
and the Press-Schechter formalism for determining the mass function, the number counts 
and the cumulative number counts as function of redshift, and study the sensitivity on the 
parameter v. We use the ACDM model to compare our predictions. Also we investigate 
how the number counts and cumulative number counts, depends on the parameters f2 m o 
and h when the value of v is fixed. 

Our paper is organised as follows. In section II we introduce the RCC model and 
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discuss the behaviour of the comoving energy density and deceleration parameter. In 
section III we discuss the linear perturbations of matter. Section IV is devoted to the 
computation of the number counts in the Press-Schechter formalism. In section V we 
present our conclusion. In the appendix we display results for the DGP model and scalar- 
tensor theory, which are useful to calculate the linear growth rate. 



2. Running Cosmological Constant 

In this section we shall introduce our cosmological model by considering a FLRW metric 
and a Universe consists of matter (baryonic + dark) and a dark energy component or 
vacuum energy; thus the cosmological evolution is governed by the following equation of 
Friedmann 

TT o 8irG , . k , 

H 2 = —(p m + PA )-—, (2.1) 



where p m and p\ are densities of matter and vacuum energies respectively, a is the scale 
factor, H = ~ is the Hubble parameter and k is a constant, which can be chosen to have 
the value +1,0,-1. In this investigation we only restrict us to the flat case. Using the 
Bianchi identity, the densities of energies have the conservation law 

dpm + dfA = 3 Pm , 2 2 , 

dz dz 1 + z 

Of the above equations we can see that there is an exchange of energy between the 
two components of the universe, so both are completely tied. However, we assume that the 
baryonic component no energy exchange with the variable CC, since that it would not be 
compatible with constraints of particle physics, but otherwise the dark matter component 
can in principle be subject to an energy exchange resulting in a time-dependent mass for the 
dark matter particles and a modification in its equation of state. Therefore, in this model, 
the dynamic of the Universe is dominated by the evolution of the CC and its interaction 
with dark matter. 



On the other hand, an important question to be noted is that the equations ( 2.1 
|2.2f ) do not form a complete set of equations, since that there are three free variables, 
Pm{ z ) i P k{ z ) i H(z) and two equations. We need a third equation to have a complete system. 
In this paper, we introduce an additional equation by considering the effects of applying 
the renormalization group to the CC. In this framework the energy density of RCC is 



22, 21 



d PA aH 2 M 2 



dlnH (4vr) 2 ' 

The above equation was proposed assuming that the renormalization group scale is iden- 
tified with the square root of the curvature scalar yi?, which in the FLRW cosmological 
context is equivalent to identify it with the H{z) at any given redshift (or cosmological 
time) . M is an effective mass parameter representing the average mass of the heavy parti- 
cles of the Grand Unified Theory (GUT) near the Planck scale, after taking into account 
their multiplicities. The coefficient a can be positive or negative, this sign dependece on 
whether bosons (<r = +1) or fermions (a = —1) dominate in the loop contributon, this is, 
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Figure 1: Comoving background matter density as a function of redshift for the RCC model. 
We can see that the decrease of density is a characteristic signature of the coupling between dark 
matter and vaccum energy. Notice that in this plot ACDM corespond to a constant line equal to 
one. From bottom to top: v — 0.25 (orange), = 0.065 (red), = 10 -3 (blue). In all model we used 
n m0 = 0.24 and h = 0.72. 



depending on whether the fermions or bosons dominate at the highest energies. Whithout 
going into the detalis on the model, let us recall that the equation ( |2.3| ) is interpreted in 
the renormalization group framework as a " /3-function" of QFT in curve-space time, which 
determines the running of the CC, for other details see ||24|]. Therefore, for background 
(cosmology of the zeroth order) these above three equations us allow a full description of 
the evolution of the RCC cosmological model. Now using the equation (|2.3| ) we can obtain 
one explicit solution for p\ 



PA = PA0 + ^M 2 (H 2 



Hi), 



(2.4) 



where pao an d Hq are the current values of these parameters and, additionally, in this 
model appears a parameter new v (dimensionless) that is given by 

a M 2 



v = 



12-k Mj '' 



(2.5) 



thus, as we can see, from the equation (2.4), if v is zero then the effects of a running are 
canceled and recover the standard model ACDM. For the case of the Planck mass, that 
is, M = Mp, the parameter v takes the value (positive or negative, depending on a), 

1 

12^ 



±2.6 x 10 



-2 



(2.6) 
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Figure 2: The deceleration parameter q(z) in the RCC model. For higher values of v faster we 
move away from a universe dominated by matter in the past. We used Q m o = 0.24 and h = 0.72. 
For v — 10 -3 (blue line) the RCC model is indistiguishable of KCDM{ red line) 

Based on the above equations the equation ( |2.1[ ) now can be written as 

* i + (»„,„ - MM qi± ')"" zl) + "«(* 2 i ^ , (2 . 7) 
Hg V i - 3^ y V i - v ) 1 - 3v ' v ; 

where #0^0(1 + z) 2 = — and f2 m0 = ^a/d + ^B0) being dark matter and f^o 

baryonic matter. This expression for the Hubble parameter is valid for Universes with 
positive or negative curvature. Therefore, the equations (2.1), (2.4) and (2.7) define our 
cosmological model. One of the first quantities that we can calculate is the matter density. 
Thus we denned the comoving matter density function as Ii(z) = 7j5p§p- in the figure 1, 
we displayed it. There is a decrease of the density due to the coupling between dark matter 
and dark energy. Increasing the coupling leads to faster decreasing of the density. In the 
figure the line constant correspond to AC DM [y = 0). Note that this behavior is similar 



to coupled quintessence models [29| and we use only positive values of the parameter v. 
To negative values the density IT is greater than one. 

Another important parameter in cosmology of background is the deceleration param- 
eter q(z) = —1 — H(l + z )-^{-jj)- I n Figure 2 we show this parameter for our RCC model. 
Here we observe that if v takes large values (positive) then it is not phenomenologically 
compatible with a matter dominated phase in the past, that is, necessary for the forma- 
tion of structures. In this sense, numerically, when v < 1CP 3 (M ~ 0.38M p ) we can see 
that the value of q(z) tends rapidly to 1/2 and is indistinguishable from AC DM. For the 
negative case q(z) is indistinguishable from AC DM for values around v > — 10~ 3 in other 
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cases the value of q(z) is greater than 1/2 in the past. In what follows we consider linear 
perturbations and structures formation. 



3. The Linear Perturbations Equations 

The theory of the cosmological perturbations is based on the expansion of Einstein's 
equations to linear order about the background metric. The first step in the analysis of 
the metric fluctuations is to classify it according to their transformations properties under 
spatial rotations. There are, scalar, vector and tensor, perturbations. In this paper, we 
only study scalar perturbations, specifically density perturbations in the RCC model, since 
that we are interested in the matter power spectrum. Other fundamental question is the 
choice of gauge or choice of a coordinate system. The General Relativity leads to the 
issue of gauge freedom. This is, metric and matter fluctuations take on different values 
in different coordinate system. A way to deal with the gauge problem is to eliminate the 
gauge dependence entirely. This approach is called of gauge-invariant variables and was 
pioneered by Bardeen [j30|. However, in the literature, many other gauges have been used 
13 31]. The more commonly used is the synchronous gauge. For example, in this gauge, the 



dynamics of density perturbations for the RCC model was investigated in [25 1 . In the 
present work we extending the analysis of by compute the matter power spectrum in 
the longitudinal gauge and to investigate the cluster number counts. 

The longitudinal gauge is from the physical point of view much more intuitive, because 
the metric perturbations are similar to the Newtonian perturbations. This gauge is com- 
monly chosen for work in CMB and gravitational lensing. On the other hand, conceptually 
this gauge fixed all freedom degree and the two scalar potentials ^ and $ that appear in 



the line element correspond to the variable of Bardeen gauge-invariant [27, |30f] , this is not 
the case of the synchronous gauge where exist a residual transformations and leads to the 
appearance of unphysical gauge. However its use is justified, since these spurious modes 
are canceled of some manner when calculating a physical observable, which by definition 
can not depend of a given system of coordinates. The metric in the longitudinal gauge is 
given by|26[| 

ds 2 = -(1 + 2$>(x,t))dt 2 + (1 + 2$(x, t))dxjdx j . (3.1) 

It should be said that the longitudinal gauge is restricted for scalar modes; nonetheless, 
it can be easily generalized to include the vector and tensor degrees of freedom [31]. Further, 



in the absence of anisotropic stress one of the Einstein's equations gives \E r = — the two 



gravitational potentials are equal and opposite [32]. Therefore, there remains only one 
free metric perturbations variable which is a generalization of the Newtonian gravitational 
potential. This justifies the name of Newtonian gauge. 

In order to derive the equations for the density perturbations we follow the stan- 



dard formalism [31, |33|, 32]. We consider the entropy perturbation as negligible and a 



energy-momentum tensor free of anisotropic stresses. Thus, in these conditions the energy- 



-6- 



momentum tensor have the form of a perfect fluid 



T$ =P9S + (p + P)U tt U v , (3.2) 

where = dx^/(— ds 2 ) 1 / 2 is the four-velocity of fluid, p is the pressure and p is the energy 
density of a perfect fluid. For a fluid moving with a small velocity v 1 = dx l /dr (peculiar 
velocity), where v l can be treated as a perturbation of the same order as bp = p — p or 
bp = p — p. The quantities p and p are with respect at the background. In linear order the 
perturbations of the energy-momentum tensor that we used are given by 

T ° = -(p + Sp), 

if = (p + p)vi = -li, (3.3) 
T) = (p + bp)b% 



the perturbated four-velocity is 



U a = ((!-*),«*), (3.4) 



while one can directly work with the Einstein's equations, it turns out to be convenient to 
use the equations of motion for the matter variables, since that, we are eventually interested 
in the matter perturbations. We consider the conservation of the energy-momentum tensor 
as 



rpOL _ . -pa rp6 , -pS rpa _ p. fn fr\ 

~ ~Q^ol <5a i a + 1 a/3 1 5 ~ U i K 6 -^) 



this expression gives us two equations a for (3 = and other for /3 = i. To determine these 
equations we consider the perturbations for the two densities and the metric 

Pm Pm (1 + 5 m ) , 

Pa^Pa(1 + S a ), (3.6) 

Indeed, in order to simplify our work we want to make full use of the symmetry under 
translations in space, which can best be exploited by working with the Fourier components 
of perturbations. Our convention, for all quantities, will be 

m = J ^j^Fik), (3.7) 



with these considerations, we use the line element (equation 3.1), the equations (|3,3j ) and 
Q3.6D into equation (|3.5|) and we get 



5pm + Sp A + p m (0 + 3tf ) + 3H8pm = 0, (3.8) 
PmO + p m 0-^- + 3H Pm 6 = -(^ftPA, (3.9) 
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where the dot is the derivative with respect to cosmic time and 6 = d{Vi. The energy 
density of the variable CC, p\, can be written as [25] 



PA = A + B(V^) 2 , (3.10) 

where we have used the fact that V ^U^ = 3H and defined that A = p^o — ^M 2 Hq and 

B = l -7^g- . Now using the perturbated four- velocity equation ( |3.4| ) , the Christoffel symbols 
perturbed and keeping only the linear order term, we find 

PA(Z) 

On the other hand, the Einstein's equations in the longitudinal gauge are 

+ 3H(z)(V + H(z)V) = -4nGa 2 (5 m p m (z) + 8 A p A (z)), (3.12) 
* + H(z)V = -4irGap m 9, (3.13) 

substituting the equation (3.13) in (3.12) we obtain an new equation for the scalar potential 
in term of variables of matter 

2 3H Pm e (5 mP + 5apa) 

2/c 2 (l + z) 3 2(1 + z) 2 v ' 

It is convenient for our numerical calculations to write the equations (3.8), (3.9), (3.11) 
and (3.14) in terms of the redshift, i.e., using dt = H ^ z +Z ^ , and also is advantageous to use 
the following two ratios |25| 



_ Pm _ (1 + z)2H{H>) - 2H 2 n kQ (l + z) 2 

Jl(z) Pt m 2 -H$n k0 (i + zy) 

PA = 3H 2 - 2H(1 + z)H' - Hfa k0 (l + zf 
MZ) Pt 3(H2-H 2 n k0 (l + z) 2 ) 



(3.15) 



and also 
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where pt(z) = 3H 2 (z) — 3i?o^fco(l + z ) 2 . Using the equations (3.15-3.17), the equations 
system for the perturbations can be written as 



Sa(z) 



vvH, 1 



3pt 



Pth'h 2P(l + z) J ' 



(3.18a) 



5 m (z) + S m (z) 
3^ 



fx , 3/i 
fi 1 + z 



+ 



+ 



1 



+ 



3p 4 



+ 



3/ 2 



(1 + *) 



5 A 



(1 + z) H(l + z) V/i fc 2 (l + z) 



(3.18b) 



7i 2 



+ 



/i 



K'(z) + 



M'(z) 



3(A - l)v k 2 (l + z) 
v (z) H r — — = ( <5a/2 



l + z 

pt 

2(l + z) 2 



if 



3f 



k 2 (l 



(3.18c) 



(3.18d) 



where the prime is derivative with respect to redshift, v = f\6 and the terms of the equation 
( 3,18| b) M{z) and K{z) are given by 

1 



M{z) 



hk^i + zy 



K(z) 



1 



(3.19) 



(3.20) 



3^/2/1 

in these equations v is the parameter defined in (2.5) and is of great importance to us, 
since when v = we can see that the perturbation on the vacuum energy is canceled. In 
this way, we recaptured the scenario ACDM as a particular case. 

3.1 The Linear Growth Rate 

The solution of the system (3.18) allows us to determine the density contrast of matter, 
5 m , which is necessary to determine the growth factor, defined as 

$m (a) 



D(a) 



6 m (a = 1) 



(3.21) 



In the figure 3 (left) we show the growth factor for the RCC model for different values 
of v and we can see that for v = 10 -4 there is concordance with a AC DAI model. How- 
ever, observationally is more important the linear growth rate that measures how rapidly 
structure is beging assembled in the Universe as a function of cosmic time (scale factor or 
redshift) and is defined by 

din D{a) 



din a 



(3.22) 



this quantity has been measured using different catalogs. In general, the redshift maps 
of galaxies are distorted by the peculiar velocities of galaxies along the line of sight. For 



-9- 



a 



Figure 3: The left figure is the growth factor and in the right figure we show the linear growth 
rate for the RCC model. In all case we use £l m o = 0.24 and h = 0.72. 



large scale such distortion can be expressed through the redshift distortion parameter /3 
and may be show to be related to the linear growth rate as 

= (3.23) 

where b^ is the linear bias value (pL = a 9 ^ 1 /u^ ass ), that is, the ratio between the root- 
mean-squared [rms) density contrasts in the galaxy and mass distributions on scale R 
where linear theory applies. Therefore, a measure of g{z) can be obtained by these two 
parameters. The f3 parameter may be measured from redshift surveys by measuring the 
power spectrum of the galaxies [34] and the &l can be obtined from the skewness induced 



in the bispectrum of a given survey. Using this technique a value of /? and bi has been 
measured using the 2dFGRS sample of 220 000 galaxies ]3q| . Recently, another measure 
was made using the spectroscopic data from the VIMOS-LT Deep Survey |J6|; and finally 
there is another measure of growth rate using the 2dF-SDSS LRG and QSO survey |}7| ; 
However, in this last case the value of /3 and bi are not fully independent, because they 
have been obtained by imposing simultaneous consistency with the clustering measured at 
z = 0. All these data are compiled in the table 1. In figure 3 (right) we show this estimates 
of the growth rate compared to predictions from various theoretical models. For plotting 
the linear growth rate for the DGP and scalar-tensor models, we use the results shown in 
the appendix. Despite the large error bars, the measurements indicate the need of a small 
value of the parameter v, and hence very close to the KCDM model. 

3.2 The Matter Power Spectrum 

Other important amount that we can determine from the set of equations (3.18) is the 
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Table 1: In this table we compile the different values of the linear growth rate that there are in 
the literature. 

Survey (3 ±5(3 bz ± 5bi g ± 5g z 



2dFGRS |35j 0.49 ± 0.09 1.04 ±0.10 0.49 ±0.14 0.15 

2dF-SDSS LRG and QSO |7j 0.45 ± 0.05 1.66 ± 0.35 0.75 ± 0.35 0.55 
VVDS 01 0.70 ±0.26 1.30 ±0.10 0.91 ± 0.36 0.80 



matter power spectrum denned as 

P(k,z) = 5 2 (k,z). (3.24) 

To set the initial conditions we shall use the BBKS approximation for the transfer 
function |38] 



T(k) = fn(l + 2.34g(fc)) 

V ' (2.34g 2 )[l + 3.89g + (16.L?) 2 + (5.46g) 3 ± (6.71g) 4 ] 1/4 ' V ' ' 

in the presence of CC, we have that q = k/(Y*hMpc~ l ), where the shape parameter is 
£ = Qmohe-^o and h = H /W0. Since observable wavenumbers are in units hMpc . We 
assume than only 4% of the cosmic density is provided by convencional baryonic matter. 
For the present time a = 1, then the matter power spectrum is usually written using the 
transfer function in the form 

P(k) = Ak n T 2 (k), (3.26) 

where n measures the slope of the primordial power spectrum (we will assume n = 1 
H ^?]]) and A is a normalization constant. In order to obtain power spectra for the our 
RCC model we evolved equations (3.18) from z = 500 to z = 0. A value of 500 ensures 
that at that time the matter is dominant. In normalizing the power spectrum followed the 
methodology presented in p5[ . 

In Figure 4 we present the matter power spectrum of our model and the power spectrum 
estimated by Percival et al., using data of the SDSS Project for Q m o = 0.24 |39[ ]. The impact 
that has the RCC model on linear power spectrum is dominated by equation (3.18a), which 
depends of the value of the parameter v. Thus, if the parameter v is zero, then, we recover 
the ACDM model, since that, 5a = 0. Consequently, a value large v entails a more 
damping of the power spectrum for the case of v positive. Meanwhile, for v negative the 
deviation of the power spectrum with respect at the observational data is even stronger. 
In principle, this feature also can be seen from the inspection of the figure 4. Theoretically 
we can understand this behavior, because in the equation (3.9) we can observe that the 
parameter v, appears as a factor of the term k 2 proportional which is derived from the 
pressure gradient. When v is negative this terms change of sign. Thus, it is expected that 
for small scales {k > 0.15) and v < 0, the pressure terms become important. In general, 
for large scales k < 0.06 (for example) the k 2 terms may be negligible. 

We can determine the mass in terms of the Planck mass for av~ 10 -5 as M~ 0.6 x 
10 _2 M p , that is, corresponding approximately to the GUT particle spectrum. Therefore, 
the RCC model is viable for smaller masses than the Planck mass. 
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Figure 4: (left) The matter power spectrum for the RCC model. The curve blue is the ACDM and 
is very close from the RCC model with the best fit v = (1.29 ± 2.39) x 1(T 5 . We used fi m0 = 0.24 
and h = 0.72. 

To get a more accurate value of v and without using the rigorous theory of parameter 
estimation, we use the x 2 statistic, that is, a quantity that qualifies the fitting of the 
observational data for a given theoretical model, and is denned as: x 2 = ( - p ° b ~ Pthe ^ , 
where and P Q b is the observational value for the power spectrum, P t he is the corresponding 
theoretical result and a Q b denotes the error bar. The index i refers to a measurement 
corresponding to given wavenumber. Thus, the probability distribution function can be 
denned as P{v) = Foe~ x2 / 2 , where Fq is a normalization constant. Minimizing x 2 we can 
determine the most probable value for the parameter v. This value is v = (1.29 ± 2.39) x 
10~ 5 , where we have assumed a fl m o = 0.24 and h = 0.72. 

Finally, in the next section we consider only v > 0, and we leave for future work a 
detailed study for values v < 0, since that in the present paper our aim is to understand 
the overall behavior of the RCC model, rather than determine observational constraints, 
which requires a rigorous application of the theory of estimation of parameters |^(| . 

4. The Number Counts 

It has long been recognized that cluster of galaxies provide a useful probe of the funda- 
mental cosmological parameters. The total abundance of cluster iV and their distribution 
in redshift should be determined by geometry of the universe and the power spectrum 
of initial density fluctuations. One of the first cosmological parameters to be constrainted 
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Figure 5: In the left figure show the evolution of the comoving volume element with redshift for 
values different of v. In the right figure we show the volume element with respect of the de-Sitter 
volume. 

was eg, the amplitude of mass density fluctuations on a scale of 8 h~ l Mpc. For example, 
recently Komatsu et al., have determined a value of as = 0.812 ± 0.026 using the data 
WMAP for the ACDM model. However, in general, the value of as is even inaccurate, 
since there is the implicit uncertainty on the value of as as a function of w (see |^T[ ) . 

Our objective in the present section is to determine the clusters number and its evolu- 
tion with the redshift and investigate whether these quantities depend significantly from v 
parameter and what value should be used for have a concordance with the ACDM model. 
For this we used the Press-Schechter (PS) p2| , f43f formalism that gave a prescription for 
estimating the mass function for a hierarchical gaussian density field. 

In the PS model the comoving number density of collapsed dark mater haloes of mass 
M in the interval dM, is given by 

dM Vvr(j(M,z)M dM FK 2(a(M,z)) 2J K ' 

where II is the comoving matter mean density of the universe and 5 C is the linearly ex- 
trapolated density threshold above which structures collapse, i.e. 5 C = 5l(z = z co \). In 
an Einstein-de Sitter (EdS) model, an overdensity region collapses with a linear contrast 
5 C = 1.686 and is the value that we adopt for our calculations. As a first approximation, 
we use this value of EdS, and we postponed for future work the use of spherical collapse 
model to calculate 5 C = Sl(z co i). In principle, this election is not far from reality because 
for homogeneus quintessence models coupled 5 C does not separate too far from the EdS 
value (see |54]]). Let us recall that the RCC model is equivalent with quintessence models 



1 23 ] . The quantity a(M, z) = D(z)aM are the rms linear fluctuations of density in spheres 



of radius R containing a mass M and D(z) is the growth factor. In our analysis the rms 
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Figure 6: The left figure shows the evolution of number counts with redshift and the effect of change 
the value of the parameter v for objects with mass within the range 10 13 < Af/(/i _1 A/©) < 10 14 . 
The right figure show the evolution of number counts for objects with mass within the range 
10 14 < M/(h~ l Mo) < 10 15 . In the both case the curve with v = 10~ D is indistinguishable from 
the KCDM (red dashed) model. We used a fixed value of f2 m o = 0.24 for all curves. 



of the smoothed overdensity is given by 

ou = (4-2) 

where M 8 = 6 x lO 14 il m /i _1 M , the mass inside a sphere of radius Rs = 8h 1 Mpc, Mq is 
the solar mass, and a$ is the mentioned above variance of the overdensity field smoothed 
on a scale of size Rs- The index 7 is a function of the mass scale and shape parameter V 



m 



1 M 

7 = (0.3T + 0.2) [2.92 + - log(— )] . (4.3) 

we use T = 0.167 [E7J and we associate galaxy clusters with dark matter haloes of the 
same mass. Our analysis of the effect of a running cosmological constant on the number of 
dark matter haloes is done by computing two quantities. The first is the number of haloes 
per unit of redshift in a given intervale of mass 

dN f ^ f Msup dn dV „, 

dfl / - — -i—j^dM , (4.4) 



dz h-K JM inj dM dzdn 

where is the comoving volume element and is given by r 2 (z)/H(z), beging r(z) = 
J Q dz/H(z). In the figure 5 (left panel) we display the comoving volume element evolution 
with redshift for the RCC model. We can see that there is a strong dependence on the 
value v. Thus, large values of v comoving volume element increases. The concordance 
KCDM model is plotted for comparison. In the right panel we plot the comoving volume 
compared to Einstein-de Sitter volume for all case. 
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Figure 7: The figure shows the integrated number counts up to redshift z for objects with mass 
M > 10 13 /i _1 M Q (left) and M > 10 14 /i -1 M©(right). For both case solid line represent v = 1(T 5 
and represents AC DM. 

The other quantity that we compute is the all sky integrated number counts above a 
given mass threshold, Mj n j, and up to redshift z 

r r°° r z dn dV 

M > M *" ] ~ L dn L, L m-Gsti** (4 - 5) 

this result is also called cumulative mass function. To compute the two quantities, we 
must choose a normalization for the number density of haloes n(M). This is commonly 
expressed by as- We choose to normalize all models by fixing the number density of haloes 
at redshift zero. This is at redshift zero all models have the same comoving background 
density II and growth factor D. Our model fiducial is ACDM (O m o = 0.24, h = 0.72) with 
a 8 = 0.9 0. 

Now we can determine the dependence of these quantities when the value v changes. 
In the figure 6, we displayed the number counts as function of the redshift, obtained from 
equation (4.4) for which there is a strong dependence on the value assumed by the parameter 
v. An increase in the value of v produces an increase in the value of Comparing the 
left panel to the right, we can see that there is a larger variation for greater values of mass. 
The results for the total number of collapsed structures above a given mass are displayed 
in the figura 7. In the left panel we plot the integration in the mass range M > 1O 13 M0 
and in the right panel for M > 1O 14 M0, in both case M sup = 10 15 . We do not use strictly 
infinite, since N(z, M > M in f) is dominated by the contribution of the lower bound of the 
mass integration range. 

The parameter v can be viewed as the coupling between dark energy (vacuum energy) 
and dark matter. Models with more coupling have higher values of dN/dz. This can be 
understood by the behavior of other observables. Seeing the expressions (4.4) and (4.5) we 
note that number counts have a dependece on the growth factor D(z), comoving energy 
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density and the comoving volume element. An increase in the comoving volume translates 
into an increase in the number counts. This effect is compensates for the decrease in 
average density and D(z). Both effects produce the observed results in the figures 6 and 
7. This behavior is similar to models of quintessence homogeneous [54|. 




Figure 8: The figures show the evolution of number counts with redshift and the Integrated number 
counts, when changing the value of the parameter Q m Oi for objects with mass within the range 
I0 13 < M/(/i _1 M ) < I0 14 (top) and for objects with mass within the range I0 14 < M/(ft _1 M ) < 
10 15 (bottom). 



4.1 Cosmological Sensitivity on f2 m o and h 

In present section, we describe how variations of the parameters f2 m o and h, affecting the 
cluster abundance. The dependence of the number counts on the cosmological parmeters 
as J7 m o and h (and also v in our case) is because the equations (4.1), (4. 4) and (4.5) depend 
on these parameters through p, the factor growth and comoving volumen. The dependence 
cosmological is implicit, but is very strong. For example: increasing in the value of the 
matter parameter directly translates into decreased of the comoving volume, this can be 
seen in the reference |4S|| . On the other hand, a decrease of the latter amount, leads to 
a decrease in cluster number counts, and of the studied quantity dN/dz that is directly 
proportional to the comoving volume. 
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Figure 9: The figures show the evolution of number counts with redshift and the integrated number 
counts, when changing the value of the parameter h, for objects with mass within the range 10 13 < 
M/(/i" 1 M ) < 10 14 (top) and for objects with mass within the range 10 14 < M / [h 1 Mq) < 10 15 
(bottom). For both case solid line represent AC DM and is indistinguishable from v = 10~ 5 . 

First we consider the effects of changing Cl m o, which are displayed in figure 8. In left 
panel are showed results for dN/dz for mass between 10 13 < M/(h Mq) < 10 14 (top 
left) and 10 14 < M/(/i" 1 M ) < 10 15 (bottom left). The curves correspond to a flat RCC 
universe with h = 0.7, v = 10~ 5 for all case, and f2 m o = 0.20 (light blue solid), Q m o = 0.24 
(dark blue dashed) and f2 m o = 0.30 (red solid). In right panel shows the total number of 
clusters N(z) for the same values of parameters. Several conclusions can be drawn from 
figure 8. Overall, for a decrease in £l m o increases the number of cluster at all redshift (and 
vice versa). The curve closest to AC DM is the central curve. Note that the dependence on 
fi m o is strong, for instance, a 16, 7% decrease in £l m o increase the total number of cluster 
N(z) by 21% for more massive structures (the bottom right panel). 

The figura 9 also demonstrates the effects of changing h. Comparing the figures 8 
and 9, the quantitative behavior of the observables (N(z) and dN/dz) under changes in 
h and Q m Q are similar: decreasing h increases the total number of clusters, but does 
not considerably change their redshift distribution for objects with mass within the range 
10 13 < M/(h- l M Q ) < 10 14 , but for mass between 10 14 < M/(h~ l M Q ) < 10 15 the change 
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is greater. 



5. Summary and Discussions 

In this paper we have determinated the cosmological implications of the RCC model using 
linear perturbations. We use the longitudinal gauge to determine the density perturbations 
of matter. In analysing these perturbations we have found, similarly as in the reference 
|25| |, that the perturbation of vacuum energy density is proportional at the v parameter 
(equation 3.18a); thus when v is zero, then is recovered the standard scenario AC DAI. 

The linear perturbations theory allows us to calculate the linear growth rate g{z) and 
compare our result with other models frequently studied in the literature and with data of 
growth rate. In left figure 3 we have seen that for values of v < 1CP 4 our model predicts 
a overdensity similar to other models (ACDM, DGP and scalar-tensor theory). In right 
figure 3 we observed that for large values of v the predicted g(z) is well off the data. These 
three data of g(z) are not sufficient to rule out a given model. In Figure 4, we compared 
the matter power spectrum with the SDSS data obtained by Percival et al. ] j39f | and showed 
that a RCC model with v ~ 10~ 5 is compatible with ACDM. The best-fit for our free 
parameter is v = (1.29 ± 2.39) x 1(T 5 . 

We investigated the expected evolution of galaxy cluster abundance in the RCC cos- 
mology. The number counts of collapsed structures are frequently proposed as a tool to 
probe dark energy models. In the present paper we have investigated, using the PS mass 
function, the modifications introduced by a running cosmological constant. Our model 
predicts a excess of sources when compared with the ACDM model. This difference is 
greater for more massive structures. The number counts predicted by our model is very 
close from ACDM for u < 10~ 5 . 

Therefore, we have shown that there is a dependence significant of clusters number 
counts in the RCC model through of the amount dark matter coupled to vacuum energy. 
Increasing the coupling between dark matter to vacuum energy, that is, increasing the value 
of v (see figures 6 and 7) increase the clusters number counts. This feature is compatible 



with the clustering properties of dark energy models of quintessence [29] . In our model, 
this effect is due to the decrease of comoving matter density and increase of comoving 
volume element when v increases. 

In general our results on power spectrum and the linear growth rate are compatible with 
a large amount of other models. For example, for the coupled quintessence model presented 



in reference |48] the interaction must be very weak in the order of 10 3 to be compatible 



with data from 2dFGRS. Recently in reference |5(]] has been studied the perturbations 
of the called model AXCDM , which includes two free parameters, the parameter v that 
has the same meaning as our case, and a equation of state parameter wx, the cosmon 



component. Both components can interact. In reference [51| was done an analysis of linear 
perturbations of these model, which includes perturbations of the dark component and the 
results were consistent with data from the 2dFGRS. 



In the reference [52] is considered a model of holographic dark energy with infrared 



decay in CDM. They use three types of cut-off (Hubble, particle and future event horizons) 
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and for all cases is displayed that there are modes of growth for the density contrast, when 
the effective equation of state is in the range —1 < w e ff < —1/3. The authors, have used a 
Newtonian approximation, and the dark component is not perturbated, being only for the 
cut-off of the Hubble horizon that the model implies a dark energy density proportional to 
the square of the Hubble parameter, similar to our case. We expect that a full relativistic 
approach of density contrast, whatever in the synchronous or longitudinal gauge, determine 
results that should deviate little from the cosmological constant. 

In reference [ 53 1 is considered a model with a cosmological term that decays linearly 
with the Hubble parameter. In that paper, the authors, consider a relativistic treatment of 
perturbations in the synchronous gauge and the dark component is also perturbed. They 
calculate the matter power spectrum and show that their results are inconsistent with 
2dFGRS data, since the matter parameter resulting is very large (fi m ~ 0.7). 

On the other hand, the references p4[ |55| explores changes in the number counts pre- 
dicted for models of homogeneous and inhomogeneous dark energy, that is, two extreme 
limits for the evolution of dark energy in the overdensity region. In the case of the homoge- 
neous dark energy density the value of overdensity inside is the same as in the background. 
In the inhomogeneous dark energy there is a collapse with dark matter. The authors show 
that there is a deviation of up to 15% with respect to AC DM in the inhomogeneous case 
and for both types of models the largest deviations are observed for massive structures 
M > 10 16 . These results are consistent with our calculations. 

Finally, in the near future, we plan to further investigate other astrophysical implica- 
tions of the RCC model. It would be interesting to consider a modification for spherical 
collapse, where the S c is function of the redshift of collapse. Also can be evaluate the 
profile of density contrast around the cluster, and supercluster of voids matter, using as 
first approximation a NFW (Navarro, Fenk and While) profile |]56| . Another important 
point is to study the concentration of haloes. For example, following the prescription given 
in reference [57] we can investigate whether the concentration of haloes in the RCC model 
decrease with increasing of the mass as in the case of the AC DM model. In principle, this 
also could provide limits on the model. 



A. Equations used to model DGP and Scalar- Tensor Theory 

In this appendix we write the equations used to determine the growth factor in the DGP 
model and the scalar-tensor theory. The growth factor is defined in equation (30) and 
obeys the equation [58, ^] 

D"(k, a) + (l + ^)D'(k, a) - J^j^ /fo a)D{k, a) = 0, (A.l) 

with the initial condition D{a) ~ a for a ~ (in the matter-dominated era). This equation 
does not take into account perturbations of dark energy and the anisotropic stress. The 
function f(k, a) expresses the connection between the metric perturbations with the matter 
density perturbations. It therefore depends on the particular gravity theory. 
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A.l DGP Model 



In the DGP model for the case of flatness and when we only have matter on the brane the 
H{z) is given by || 

H 2 DGP = + \AWi" 3 + a- c , (A.2) 



where Q rc = — H m o) 2 and for this theory the f(k,a) function is given by [ |5S| ] 

f(k,a)= + (A.3) 



with 



H Qy /VL rc 3H DGP (a) 



A.2 Scalar- Tensor Theory 

Scalar-tensor theory is simplest generalization of the General Relativity where the funda- 
mental constants can be variable. The function f(k,a) in this case is given by [58| 

fl ) = r Ge i [a \^ ( 1 + TTTk~i — s) ' ^ 

G e ff(a = l)\ \ + {k/ma)J 

where G e ff is the effective Newton's constant when the scalar factor is a and a = 1 is the 
present value, while m is the mass of the scalar field $ inducing a Yukawa cut-off to the 
gravitational field. We used a simple ansatz [59] 

Geff (a) 



Geff (a = 1) 



! + £(! -af. (A.6) 



In the case of the DGP model and scalar-tensor theory the detection of an f(k, a) ^ 1 
would be a signature of alternative theories of gravity. For the numerical calculations we 
used £ = -0.2. 
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